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ABSTRACT 

We consider some applications of the Renormalization Group flow equations ob- 
tained by resorting to a specific class of proper time regulators. Within this class 
a particular limit that corresponds to a sharpening of the effective width of the 
regulator is investigated and a procedure to analytically implement this limit on 
the flow equations is shown. We focus on the critical exponents determination for 
the O(N) symmetric scalar theory in three dimensions. The large N limit and 
some perturbative features in four dimensions are also analysed. In all problems 
examined the results are optimized when the mentioned limit of the proper time 
regulator is taken. 
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1 Introduction 



In the past years there has been growing interest for the Renormalization Group 
(RG) techniques originally inspired to the Kadanoff- Wilson blocking procedure 
H , as they represent promising means to investigate areas that are out of reach of 
the standard perturbative quantum field theory. The central idea is the descrip- 
tion of the RG flow in terms of differential equations in the momentum space. 
Some of the original works on this subject are |2], f|, More recently, a new 
formulation of the problem, usually called Exact Renormalization Group (ERG), 
has been developed || [7|, H (for recent reviews on the subject and on various 
applications see||, |10[)- The ERG equations are functional differential equations 
which describe the momentum scale k dependence of a particular functional, usu- 
ally indicated as average action. This functional is in some sense similar to the 
effective action, which is the Euclidean generator of the one particle irreducible 
correlation functions, with the difference that in the former only fluctuations 
with momenta larger than k are included. This separation between high and low 
frequencies is realized by the introduction of a smooth cutoff regulator R k . A 
relevant property of this functional is that it interpolates between the classical 
Euclidean action and the full effective action of the theory, when the scale k 
is lowered from the ultraviolet cutoff, where the theory is defined by the local 
classical action, down to zero where all the fluctuations are integrated. 



A slightly different approach to the RG flow has been formulated in | TT . IT2L IE 
where an operator cutoff is introduced by means of the Schwinger proper time 
regulator. The most remarkable feature of this regulator is that it is formulated in 
a gauge invariant way and it has been already employed to compute the one loop 
beta function of Yang-Mills theories ||14j|. Therefore it is a promising tool for any 
analysis involving gauge theories. Besides that, it has been used to derive flow 
equations for scalars coupled to fermions||15|| and for O(N) scalar theories [|16|j . 



The differential equation that describes the flow of the scale dependent action 
obtained via the proper time regulator is 

dT k 1 roo ds df k , 5 2 r M 

with the constraints on the dimensionless cutoff function fk{x) that /fc(0) = 1 
and fk(x — > oo) = and the initial condition that at some ultraviolet scale A 
Tfc is equal to the classical action S. The trace sums over all the discrete and 
continuous indices of the field <p. A connection to the perturbative one loop 
effective action 

I— = -iTr f ± exp (sfpi (2) 
2 Jo s ^ ' " 
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is obtained if one neglects the effects of the running by freezing in the r.h.s. of 
Eq. (HD at its initial condition S. In this case T k , due to the condition /^(O) = 1, 
coincides at k = with Y 1 ~ loop . The running quantity Tk in Eq. (|1|) corresponds 
to an infinite resummation of diagrams and therefore improves on the one loop 
approximation. In practice the regulator kdkfk acts as a particular smooth cutoff 
or weight function in the integration of the modes, selecting only the modes within 
a small shell centered around the scale k. 

However there is no direct connection to the smooth cutoff Rk, introduced in 
the ERG formalism (see for instance H). As a consequence, the proof that in 
the limit k — > the ERG flow of the average action Tk converges to the usual 
effective action cannot be simply translated to the flow determined by the proper 
time regulator and a definite indication that this flow actually converges to the 
full effective action is missing. 

It should not be forgotten, however, that for practical purposes it is impossible 
to deal with the flow of the full action and, typically, a semi-local derivative ex- 
pansion of the action is introduced to reduce the problem of the evaluation of the 
flow to a treatable set of coupled partial differential equations. The best approx- 
imation analysed so far is the next to leading order in the derivative expansion 
which corresponds to take the following ansatz for Tk 

r fc [0] = / d D x{v{<j>) + ( 3 ) 

neglecting all the higher derivative terms, and determine the corresponding differ- 
ential flow equation for the potential V and for the wave function renormalization 
Z (obviously V and Z depend on the scale k but this dependence is not explicitly 
indicated to simplify the notation). Due to this approximation, some uncertainty 
is introduced about the k — > limit of the truncated Tk and a scheme dependence, 
basically related to the specific regulator employed, appears in the determination 
of the various physical quantities. This unphysical feature, which would disap- 
pear if one could deal with the flow of the full action, cannot be easily quantified. 
According to this remark, it is worthwhile to compare the predictions obtained 
with the proper time regulator (although lacking of a convergence proof) with the 
results obtained from the ERG. In |T7| it was noticed that the following choice of 
the function fk 

i=Q % - 

parametrized in terms of the integer m, when the value of m is sufficiently large 
provides a determination of the critical exponents rj, u, to for the one-field scalar 
theory in D = 3 dimensions that is certainly comparable to the other determina- 
tions obtained in the ERG framework. 
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In this paper we shall follow this point of view and explore the possibility of 
explicitly studying the limit in which m is infinite. It will be shown that for a 
simple particular parametrization it is possible to take the limit m —>■ oo and get 
two coupled equations for the potential and the wave function renormalization 
which do not contain the parameter m anymore. As noticed in |17| for growing m 
the proper time regulator selects smaller and smaller momentum shells, centered 
around the scale k, which provide a relevant contribution to the flow differential 
equations. In this sense the limit m — > oo corresponds to a procedure of reducing 
the characteristic width of a smooth cutoff and taking the sharp cutoff limit. Due 
to the singular nature of this limit, the operation is not unique and depends on 
the specific smooth cutoff chosen. The central result of this paper is that the 
equations obtained when m — > oo yield a determination of the critical exponents 
which actually agrees very well with the extrapolation of the corresponding values 
obtained at finite m, thus confirming that these universal quantities for large m 
do not trivially diverge or vanish but rather have a finite sensible limit. Moreover 
these values are very well compatible with other determinations stemming from 



different, well established techniques (see |18[ or In Sect. 2 we consider 

the coupled flow equations of V and Z for a scalar theory obtained with the 
proper time regulator and determine their form in the limit m — > oo. As a simple 
application we show that the perturbative two loop anomalous dimension of the 
scalar field is easily recovered. In Sect. 3 we collect the numerical results on the 
critical exponents at the non-trivial fixed point in D = 3 dimensions for finite 
values of the parameter m and in the limit m — > oo. In Sect. 4 the particular 
case of the large N limit of the potential equation for a scalar O(N) symmetric 
theory is considered in order to check that the proper time regulator reproduces 
the known exact results. The conclusions are summarized in Sect. 5. 



2 One component scalar theory 

Starting from the evolution equation (fj), with the ansatz (Q) for the action of a 
single scalar field theory in D dimensions one can extract the flow equations for 
V and Z. The procedure is outlined in [[17] and the result is 

, dv up I Zk2 Y +1 ' D/2 n 

k M =ak [zk^ (5) 
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dz D zk 2 

k—— = ak ' 



m+l—D/2r, 1 n /o \ 

Zfc 2 + V 7 ' V 



(4 + ISD - D 2 ){Z'f\ (10 - D)(m + 1 - D/2){m + 2 - D/2) 



24Z / 6(Z/c 2 + V"") 2 



(m + 1 - D/2)(m + 2 - D/2){m + 3 - D/2) ^, T//// x 2 



6(ZA: 2 + V") 



(6) 



where each prime indicates a derivative w.r.t. the field and the constant a is 
expressed in terms of gamma functions T(x) 

r(m + 1 - D/2) 

(4vr)^/ 2 r(m + 1) 1 J 

One can immediately realize that the quantity V" / {Zk 2 ) plays a fundamental role 
when the large m limit of Eqs. (|5|,|6|) is considered. In fact if V" / {Zk 2 ) > the term 
raised to the m-th power in vanishes, and if V"/ {Zk 2 ) < it diverges (we 
are not interested here in the broken phase sector and V" /{Zk 2 ) < —1 is never 
considered). This however does not imply that the universal properties related 
to Eqs. (|5|,|B|) must have a singular (or trivial) behavior for large m. Actually in 



17] it is found that already for values of m around 10, the critical exponents are 



almost m- independent. 

An interesting information about V" can be read in Fig.l where the particular 
case of the non-Gaussian fixed point in D = 3 for the potential only (and Z = 1 
fixed) is examined (more details about this calculation are given in Sect. 3). The 
triangles in the upper frame of the figure correspond to the values of —V"{(f) = 0) 
at the fixed point, properly expressed as a dimensionless variable (see Sect. 3), 
plotted for some values of 1/m (note that in this particular case the curvature at 
the origin is negative). We consider 1/m instead of m because we are essentially 
interested in the large m limit. The circles in the upper frame correspond to the 
product —mV"(<fi = 0). In the lower frame the circles are the same of the upper 
part of the figure, but plotted on a magnified scale. It is evident from Fig. 1 that 
V"{cf) = 0) at the fixed point behaves, at least in first approximation, like 1/m 
for large m. Of course this does not simply mean that mV"{(j) = 0) is a constant. 
Rather, the plot in the lower frame suggests a behavior like —mV"{(j) = 0) = 
a — b/m with a,b positive constants (up to the accuracy considered). However, 
our aim here is just to find out the leading behavior of V"(<f> = 0) for large m, 
which, as noticed above, is 1/m. Supposing that V" has the trend of Fig.l, not 
only at <fi = but at any value of <p, it would be helpful to redefine the variables 
in Eqs.(|5],|6p in such a way that this dependence on m could be compensated. 
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As a first step we can get rid of the constant a in Eqs. (f^lf) by simply 
redefining V — > aF, — > y/a<j) and leaving Z unchanged. According to Eq.(|3]), 
this operation amounts to rescale the action by the constant a. Then, following 
the indications of Fig.l, we can make a second transformation by replacing <fi — ► 
y/m<p. This time the action is not uniformly modified. With the help of the latter 
substitution we can explicitly take the m — > oo limit of Eqs. (^]||) and get 

k ^ = k D e -V"/(Z^) (g) 



k— = ],D e -V"/(Zk 2 ) 

dk 

( Z" (4 + 18D - D 2 )(Z') 2 (10 -D)Z'V" Z(V 



m\2 " 



Zk 2 24Z 2 k 2 Q(Zk 2 ) 2 Q(Zk 2 )\ 



(9) 



We see that the terms raised to the m-th power in (|5],|6]) do not turn into singu- 
larities, but instead into exponential terms. In addition, other terms of Eq. flf]) 
do not appear in Eq.(|S|) because they are suppressed for m — >• oo. Eqs. (|8],|S]) can 
now be used to study the critical properties. 

A simple check on Eqs. (^J^) and (|8|,|9|) can be performed by looking at the 
perturbative determination in D = 4 of the anomalous dimension 77 of the field 
0, defined as 77 = —kdkln(Z) where Z is the residue of the mass pole of the 
two point Green function of the theory. The specific truncation of the derivative 
expansion of the action above considered is actually sufficient to derive the lowest 



order non-vanishing perturbative contribution to 77 (see for instance [20| where 
the same calculation, in a different framework, namely starting from the Wegner- 
Houghton flow equation, has already been performed). In fact if one starts with 
the classical Euclidean action (we require the invariance under the transformation 

0^-0) 

s= 1 -(d<p) 2 + ^- (10) 

then the lowest order values, in the perturbative parameter A, of Z and V in 
d^,|), are respectively 1 and A0 4 /4!. If these lowest order values are inserted in 
Eq. (0) and the function Z is conveniently expanded in even powers of the field 



z 2 (p 2 Z 4 (j) 4 

Z + — c T- + — 77- + ■■■ (11 



then, by comparing the terms proportional to <fi 2 in (^), one immediately sees 
that Z2 is 0(X 2 ). More precisely, by integrating the differential equation from the 
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ultraviolet cutoff A to the scale k << A (but still k much larger than any other 
infrared scale such as the field mass pole) 



Consequently, Zq, which at the leading order is z = 1 (as it follows from (|10|) 
and (|IID ), to the next to leading order gets a 0(A 2 ) contribution from z 2 that 
enters the flow equation of zq through the term Z" in Eq. (Rh. We get 
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The two point Green function is obtained from the double functional derivative 
of (H) at (j) — 0. It is easy to realize that to the next to leading order in the 
derivative expansion we get T^(p 2 ) = zop 2 + V"($ = 0) and therefore, to this 
order, Z = z . It follows 

i = ^(ir*V + °<*'> < 14) 

The important result here is that the perturbative value of rj (see e.g. fl"8|j ) 
corresponds to the large m limit of Eq. (14), which is in agreement with our 



conjecture that the physically relevant results are obtained for m — » oo. 

One remark is to be made here about other terms neglected in the derivative 
expansion (|3]) that are 0(A 2 ). In fact terms like Y <p 2 d n (f)d n (j) with n > 2, which are 
neglected to order 0(d 2 ), can be computed with the same method employed to 
determine Eqs.(^|,pJ)and it turns out that Y oc X 2 /k 2n . However the contribution 
of Y to r( 2 )(p 2 ) is proportional to \ 2 {p 2 /k 2 ) n and, as long as we are interested 
in a momentum range where k 2 is much larger than a mass pole, then k 2 » p 2 
and, consistently, we can neglect the contributions of Y and similar terms which 
are suppressed by the factor (p 2 /k 2 ) n . 

Finally we turn to Eqs. (^,0) and perform the same perturbative analysis. In 
order to get the same normalization used to derive Eq. (|14"D , one has to remember 
that, globally, the following substitutions have been performed to go from (|5||J) 
to (|8|,|9|): V — > aV, — > y/am(p and no rescaling on Z. As a consequence, instead 
of starting with the classical action S in fliCf ) we have to start with a modified 
action S, equal in form to (|10D , but with A replaced with A = am 2 A, which for 
m — » oo becomes (see Eq. (^)) 

1D7T 2 
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Note that since Z is unchanged the initial value in flTT^), Z = 1, i.e. z Q = 1, 
is unchanged too. A non-vanishing initial value for z 2 would have required a 
rescaling according to the definition ([11]). After this remark, it is straightforward 
to derive z 2 and rj (through z®) from Eq. (^|): 

z 2 = ^ 2 +0(X 3 ) + 0(l/A 2 ) (16) 

V = j + 0(X 3 ) (17) 

As expected the large m limit of Eq. ([14]) coincides with Eq. ( ]T7|) which is the 
perturbative value of rj. 



3 Numerical results 

In this Section we focus on the problem of determining, through a numerical 
analysis, the fixed point solutions of our flow equations and the critical exponents 
v and uj related to the eigenvalues of the linearized version of the flow equation 
around the fixed point. In particular we are interested in the three dimensional 
case for which a large number of determinations of these universal quantities, 
obtained employing many different techniques, is available in the literature (see 
e.g. (ffSl [19], §]). Therefore, for the rest of this Section we set D = 3. 



The numerical determination of the exponents for this particular kind of 
proper time regulated flow, was started in |17| |. There, z/, oj and the anomalous 



dimension rj of the one component scalar theory, determined by the numerical 
resolution of Eqs. (§) and (||) are displayed for increasing values of the parameter 
m and the data show a monotonous dependence on m. In particular v and uj 
are determined up to m = 9 and rj up to m = 40. Since in the second case m is 
sufficiently large, the asymptotic value (for m — > oo) of r] is obtained through a 



numerical fit in [T7[ and it is only claimed that v and u converge to finite values. 



In the following we refine this analysis and evaluate 77, v and uj to the leading 
order (0(d°)) in the derivative expansion (by solving the potential part of the 
flow keeping Z = 1 and 77 = fixed), and to the next to leading order (0(d 2 )) 
(by solving the full coupled problem for the potential and the wave function 
renormalization) . In most cases we push the parameter m up to m = 120, which 
clearly shows the numerical trend of the exponents and then we compare the 
obtained values with the ones determined from the equations (§|) and ([|) which 
hold for m — > 00. We also extend the analysis to the O(N) symmetric theory 
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limiting ourselves to the 0(d°) approximation in the derivative expansion and 
derive v and u for iV = 2, 3, 4, for increasing values of m and for m — > oo. 

In order to perform the numerical analysis of the flow equations to find the 
fixed point solution and the related critical exponents we need to reformulate Eqs. 
(H) and (^) in terms of dimensionless variables. At the same time, to simplify 
the numerical procedure, it is convenient to get rid of the constant a through 
the transformation already used in the previous Section and, more important, to 
replace the differential equation for the potential with the corresponding equation 
for the potential derivative w.r.t. the field. In fact, as already noticed in || Q 
the former equation is stiff and therefore more difficult to treat. To this aim 
we introduce the new dimensionless variables t = ln(A/fc), x = k~^ +ri '' 2 (j)/\/a, 
Vd(t,x) = k~ 3 V(k,4>)/a and Z d (t,x) = k v Z(k,(p) where the subscript d has been 
used to indicate the dimensionless potential and wave function renormalization, 
a is defined in Eq.(^) and, as before, r\ indicates the anomalous dimension of the 
field. Finally we define the potential derivative f(x,t) = d x Vd(t,x). In terms of 
these new variables Eqs.([5]|)D read 



9tf JA^l f Jlpl xf 



m 




m-3/2 



z d 'f - z d f") 

(z d + n 



d t Z d 



-vZd-^xZd^ 



7(m 2 - 1/4) 
6(Z d + f') 2 



Zd'f" + 



m-l/2 r 



fm-1/2)/ „ A9(Z d '" 2 



\Zd + f'J 
(m 2 - 1/4) (m + 3/2) 

Q(Z d + f'f 



z d (f 



n\2 



2AZd 



(19) 



Analogously, we can rewrite Eqs. (^) and ([|) in the dimensionless form 

a ^(l^ f -(l^ xf + e -^(il-fj!f} (20) 



d t Z d = -rjZ d - 



(i + v) 



xz d ' + e-r'« - _l^ + im (21) 

d+ [ Z d 24Z 2 6Zj + QZj ) {21) 



The fixed points are the t- independent solutions f*(x) and Zd*{x) of Eqs. 
18|jl9]) , or Eqs. ( ]20| , |2T|) . The resolution of the coupled ordinary differential fixed 



point equations also provide the value of r\ at the fixed point. In order to find 
v and u we need to introduce small t-dependent perturbations around the fixed 



S 



point and reduce the problem to coupled linear equations in the perturbations. 
The perturbations are parametrized in the following way 



V i (t,x) = V d *(x) + e xt v(x) 
Z d {t J x) = Z/(x) + e xt z(x) 



(22) 



The corresponding relation for the potential derivative f(t,x) is obtained by 
deriving the first line and, by indicating with h(x) the derivative of v(x), one has 



f(t,x) = r(x)+e M h(x) 



(23) 



The resolution of the linear equations provides the eigenvalues A and the eigen- 
vectors v(x), h(x) and z(x) which govern the flow close to the fixed point. In 
particular for the problem considered we expect, besides the Gaussian fixed point, 
which corresponds to f*(x) = n = 0, Z d [x) = const, only the Wilson-Fisher 
fixed point, which has just one relevant eigenvector (A positive). The exponent v 
is defined as the inverse of the only positive eigenvalue and u is the opposite of 
the less negative eigenvalue. 

We shall also consider the O(N) theory but only in the lowest order approx- 
imation 0(d°), and in this case the dimensionless form of the flow equation for 
the derivative of the potential is 



dtf — —f x f ' + (m 

tj 2 j 2 V 2, 



+{N 




m+l/2 



1 + f/x 



f 



x- 



(24) 



where in addition to the longitudinal field contribution, which is equal to the one 
for the single component field, there are (N — 1) contributions of the transverse 
fields. Finally the limit m — > oo of Eq. ([24]) becomes 



dtf = -f 



\xf' + e-rr + {N-l)e-"* 



(25) 



It must be noted that Eq. ( p4|) at m = 1 and for any iV is exactly equal to 
the flow equations derived in the 0(d°) approximation in |8], 21] and therefore a 
good check on our numerical outputs, at least for m — 1, is to compare them to 



the results reported in [21 



The numerical resolution of the equations is thoroughly explained in || and 
here we just shortly summarize the procedure and skip the details. The two 
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coupled ordinary differential equations that determine the fixed point have some 
boundary conditions defined at x = and the rest at large x. In fact, by re- 
quiring the symmetry of the theory under the transformation x — > —x, one gets 
the following conditions: Z2'{0) = f*(0) = 0. In addition, the normalization 
Z^(0) = 1 is imposed. (As noticed in the flow equations here considered are 
reparametrization || [22[ invariant, and, as a consequence the universal quantities 
are independent of the particular normalization chosen, and this has been checked 
explicitly in the numerical analysis of ]17|])- On the other hand, the equations 
themselves constrain the asymptotic behavior of the solution at large x, up to 
two constant factors, one for each of the two functions Z%(x) and f*(x). The 
eigenvalue problem for the linear equations is totally analogous because the same 
symmetry must be imposed at x = and the same kind of constraints can be 
derived from the equations at large x. 

These boundary conditions are sufficient to solve the problem (and determine 
the fixed point and the critical exponents) by means, for instance, of the shooting 
met hod . In particular, for our analysis, we have implemented this method by 
making use of the NAG libraries and this has allowed us to improve on the 



accuracy of the results of [1 7 and to test the equations for higher values of m. 



One result of the numerical analysis has already been presented in Fig. 1, 
where the curvature at x = of the non-Gaussian fixed point solution for some 
values of m, is plotted. Namely the triangles in the upper frame of Fig. 1 are 
the values of —10 f*'(0) determined from Eq. (|I8| ), in the 0(d°) approximation 
i.e. with Zd = 1 and rj = 0, plotted versus 1/m. The circles in the upper and 
lower frames correspond instead to the product —10 m/*'(0) displayed with two 
different scales on the vertical axis. This plot has already been commented in 
Sect. 2. 

Fig. 2 shows the anomalous dimension rj at the fixed point for the N — 1 case, 
plotted versus 1/m up to m = 120. Note that we have inserted at 1/m = the 
value of rj obtained from the fixed point solution of Eqs. (^,^1|) whereas all the 



other points come from Eqs. (|18yi9l) . Some of these points are already reported 
in [0, where however the maximum value of m considered was m = 40. The 
convergence in Fig.2 of the values obtained at finite m to the one at 1/m = is 
clear. 

In Fig. 3 we have collected the same kind of results for v and u. Namely, for 
the N = 1 case, empty and black circles are respectively the O(d ) and 0(d 2 ) 
values of v and empty and black triangles the same for u. Again the points on 



the vertical axis 1/m = are obtained from Eqs. (pn,Bl|) and the others from Eqs 



(|T8|JT9|) . Note that for the O(d ) case we evaluated the exponents up to m = 120 



whereas for the 0(d 2 ) we stopped at m = 10 because we believe that even at this 
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level the trend of the data and their convergence to the 1/m = case is evident. 



Two remarks are in order. The first one concerns the fact that, for those 
values of m for which the comparison is possible, we have found agreement with 
the numerical data of ]T7|, except for u to the 0(d 2 ) order, whose values plotted 
in Fig. 3 (black triangles) differ from those of In fact, we have noticed that 
in this case the results are more sensible to the maximum value of x considered in 
the numerical integration of the equations and the routines used here allowed a 
more careful analysis and a more accurate determination of the data. The second 
remark is about the convergence of the derivative expansion. The 0(d°) and the 
0(d 2 ) values of v at 1/m — are almost equal, suggesting that the derivative 
expansion converges very rapidly for this particular exponent. The situation for 
u} is rather different, and one can only guess that the higher order determinations 
of this quantity lie between the 0(d°) and the 0(d 2 ) estimates. Therefore it 
seems that the rapid convergence concerns only the relevant eigenvalue. 

The findings of the numerical analysis for v and uo in the N = 2 theory and 
0(d°) approximation are also reported in Fig. 3. The symbol + is associated to 
the values of v and the x to uo. In this case we have determined the exponents 
up to m = 120 and again the convergence toward the values obtained at 1/m = 
is clear. On the basis of these results, we do not consider so many values of m 
for iV = 3, 4, and determine v and to only at m = 1,2 and at 1/m = (see the 
Tables below). 

In Fig. 4 we report various plots obtained from Eq.(^), which is valid to 
the 0(d°) order, in the four cases iV = 1,2,3,4. In the upper frame, the fixed 
point derivative of the potential is shown, in the central frame the eigenvector 
corresponding to the critical exponent v and, in the lower one, the eigenvector 
corresponding to uj. The trend observed is a flattening of the curves for increasing 
N. 

Finally in Table I and II we report the numerical values of the exponents only 
at m = 1, 2, oo. Incidentally we note that rj at m = oo in Table I is practically 
the same of the value rj = 0.0329 quoted in [T7]] as the output of a fit to the data 
collected up to m = 40. We also note that the 0(9°) results at m = 1 are totally 
in agreement with those of GTJ which, as discussed above, was expected because 



in this particular case the flow equations considered here and in |21| coincide. 



A comparison with the critical exponents determined by different methods (see 
19, ||) globally indicates a better agreement with the values at m = oo rather 



than with those at m = 1 or m = 2. This supports the idea that the optimization 
of the proper time cutoff is achieved in the large m limit, which can be interpreted 
as a sharp limit of this class of regulators. Conversely, the optimization criterion 
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suggested in P5|, obtained by establishing a relation between the proper time 



regulator and the Rf, class of regulators introduced in the ERG approach, does 
not seem to be particularly satisfying. In fact it predicts that the optimum value 
of our parameter (note the difference in the definition of m here and in p3[ ) is 
m = 3/2 and the exponents at m = 3/2 are expected to be between those found 
at m = 1 and those at m = 2 in Tables I, II. 



4 O(N) theory at N = oo 

Another interesting check on the proper time flow equations concerns the problem 
of the O(N) symmetric theory in the limit iV — > oo. In this case the O(N) theory 
reduces to the spherical model which is exactly solved. At the same time in this 
limit the flow equation for the potential becomes a close decoupled equation and 



the main universal features of the spherical model are recovered |2], |24|, |2T|. Our 
aim is to check that these properties survive when we interchange the two limit 
procedures and consider first m — ► oo and then iV — > oo. As noticed before, 
to the lowest order in the derivative expansion, the potential flow equation (^4j) 
at m — 1 is equal to the one considered in |2~I|, where a detailed analysis of 



the N = oo case is presented. So at least for m = 1 the proper time flow of 
the potential reproduces the results shown in ETJ and in particular the presence, 



besides the Gaussian fixed point, of the Wilson-Fisher fixed point for D < 4, and 
the correct determination of the eigenvalue spectrum of the linearized equation 
in both cases. Then we examine the problem at m ^ 1. 



Let us consider the more general version in D dimensions of Eq.(23), which 
was obtained in D = 3. For large N, the leading behavior is given by the last 
term in the r.h.s.. As it has been done for the constant a, we can make one more 
rescaling / — > Nf and x — > y/~N x which allows to take straightforwardly the 
limit N —>■ oo 

df ( D + 2 m )f ( D ~ 2 m ) xf > (26) 

tJ \ 2 x 2 {l + f/x) m+1 ) J \ 2 x 2 {l + f/x) m+1 ) J y ' 

The fixed point structure is obtained by requiring dtf = 0. Actually we are 
interested in the universal properties at the fixed point and, as shown in [0, p4]j , 
to determine the critical exponents it is not necessary to have the explicit fixed 
point solution f*(x), but it is sufficient to determine the essential singularities of 
the equation. From the structure of Eq. (p6|) one realizes that there is a non-zero 
value of x at which the function /* vanishes. In fact if at x = x the second bracket 
in (|2~6D vanishes, then it must be f*(x) = since the first bracket at x = x does 
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not vanish. Therefore from the second bracket it follows that x = y2m/ (D — 2). 
We can now expand around x. By setting in the fixed point equation x = x + Ax 
and f*(x) = f*'(x)Ax + 0((Ax) 2 ) we get the value of f*'(x): 

( m + !)/*>) = _±_- 2 (27) 

With this information we can analyse the corresponding linearized equation and, 
recalling Eq. (|2~3"D, we get 




ni 



x 2 (l + f*/x) m+1 
m 



h 



-) X h' + m{m+1){r - Xr,) h (28) 
x 2 {l + f*/x) m+1 ) x 3 (l + f*/x) m + 2 y ' 

At x = x the coefficient of xh' in brackets vanishes and therefore from Eq. 
(|28| ) it follows that either h'(x)/h(x) is singular (with a simple pole singularity) 
or it is finite and in this case the eigenvalue A is uniquely determined (remember 
that f*(x) = 0). If one requires that the eigenvector h(x) has no singularity and 
admits the power expansion around x, 

oo 

h(x)=^£ i a i (x-x) i (29) 

i=n 

where the lowest power n is a non-negative integer, then one can study the behav- 
ior of Eq . (|28|) close to x = x. The quantity h'(x)/h(x) is replaced by the leading 
singular term n/(x — x) but the pole is canceled by the vanishing coefficient of 
h'(x)/h(x), and therefore with the help Eq. (|27|), at x = x Eq.(p8|) is reduced to 
the simple relation among the eigenvalue A, the non-negative integer n and the 
number of dimensions D 

A = D - 2 - In (30) 
which is the spectrum of the spherical model eigenvalues. 

We are also able to derive the eigenvalues relative to the Gaussian fixed point. 
In fact, going back to the linearized equation (|28|) , by evaluating it at the Gaussian 
fixed point fc(x) = 0, and expanding around x, one gets the new spectrum 

\ G = 2-(D-2)n (31) 



At the upper critical dimension D = 4 the two spectra (^) and (|3lD coincide. 
Obviously for the simple case of the Gaussian fixed point one could directly start 
with the linearized equation of the potential, rather than using the potential 
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derivative /. By indicating with v^ the eigenvector of the potential, the N = oo 
flow equation for the potential, linearized around the Gaussian fixed point Vq = 0, 
rccicls 

(A - D), c = g - <^z) v c < (32) 
The above equation is solved explicitly 



2m 



(D-X)/(D-2) 



VG = c [ x2 -W^Y)) (33) 

where c is the integration constant. If we require polynomial (not constant) 
solutions of the problem, which preserve the symmetry x — > —x, we have (D — 
Xg)/{D — 2) = n + 1 where n is again a non-negative integer. Therefore we 
have found again relation (|3~TD. It should be noted that here we have the full 



solution (|33"D of Eq. (j3~^ ) and its parity must be explicitly required. Conversely 
this constraint must not be applied to the expansion considered above, around 
the point i/O. 

The relevant point of the above analysis is that in both eigenvalue spectra 
for the Gaussian and the non-trivial fixed point, the dependence on the proper 
time regulator parameter m has disappeared and the structure of the eigenvalue 
spectrum is valid for any value of m. We expect that by performing first the m — > 
oo and then the N — > oo limit in the flow equation, these universal properties are 
preserved. To prove this point we perform on Eq. (p5|) the same rescaling used 
to derive Eq. (EH) and get 



, /H ^)/-(^VW^] (34, 



As before, in the fixed point equation there is a value of the dimensionless field, 
x = \J2/(D — 2), such that f*(x) = and, as can be seen by expanding the fixed 
point equation around x, f*'(x) = 2(4 — D)/(D — 2). Then we can analyse the 
linearized version of Eq.(|34|) around the Gaussian and the non-Gaussian fixed 
point as we did for Eq.(^) and it is straightforward to derive again the two 
equations for \q and A, ( |3"T| ) and (§(]), respectively. Therefore, even in this case, 
the flow equation obtained in the limit m — > oo retain the relevant universal 
properties observed at m — 1. 
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5 Conclusions 



As already discussed in the Introduction, there is no proof of convergence of 
the proper time RG flow to the full effective action when the infrared scale k 
is lowered to zero, so that all modes in the momentum space are effectively 
taken into account. On the other hand we have collected here some evidences 
showing the reliability of the predictions of this flow equations which are certainly 
comparable to those of the ERG. 

The proper time regulator is parametrized by a number m, which is restricted 



to integer values. According to [17], m is related to the effective width associated 
to the proper time cutoff and larger values of m correspond to more narrow 
widths. We have shown that for a particular redefinition of the quantities entering 
the flow equations, it is possible to formally take the limit m —>■ oo and get new 
differential equations which are totally m independent. 

Then, in some numerical examples, we checked that the values of the critical 
exponents rj, v and u as functions of the parameter m are effectively converging 
to the numerical values obtained from the asymptotic (in the sense of m —>■ oo) 
equations. As a consequence of this analysis we get an optimized determination 
of the above exponents for m — > oo which are certainly in very good agreement 
with the numbers coming from different methods. Other aspects investigated are 
the perturbative determination of the anomalous dimension in D = 4 and the 
comparison with the spherical model in the large N limit. In both cases the 
asymptotic equations produce results that totally agree with the ones coming 
from the m dependent flow equations at large m. In the N = oo case, actually, 
the universal quantities show no dependence at all on m. 
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TABLE I 



N= 1 




7] 


V 


uj 


m 


o{d 2 ) 


0(9°) 


0[d 2 ) 


0(9°) 


0[d 2 ) 


1 


0.0653 


0.6604 


0.6337 


0.629 


0.677 


2 


0.0507 


0.6439 


0.6304 


0.674 


0.700 


oo 


0.0330 


0.6260 


0.6244 


0.762 


0.851 



Table 1: Critical exponents for N = 1 at the Wilson-Fisher fixed point to the 
leading, 0(cP), and next to leading order, 0(d 2 ), in the derivative expansion. 



TABLE II 





N = 


= 2 


N = 


= 3 


N = 


= 4 


m 


V 


uj 


V 


UJ 


V 


UJ 


1 


0.7252 


0.662 


0.7809 


0.707 


0.8238 


0.751 


2 


0.6990 


0.681 


0.7498 


0.700 


0.7924 


0.726 


oo 


0.6689 


0.746 


0.7102 


0.734 


0.7481 


0.729 



Table 2: Critical exponents v and uj for N = 2,3,4 at the Wilson-Fisher fixed 
point to the leading order 0(d°) in the derivative expansion. 
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Figure 1: Upper frame: —10 /*'(0) (triangles) and —10 mf*'(0) (circles) at 
various values of 1/m, for iV = 1 and to the leading order in the derivative 
expansion. Lower frame: The same circles of the upper frame, plotted with a 
magnified scale in the vertical axis. 
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Figure 2: Anomalous dimension rj at various values of 1/m, including 1/m = 0, 
for N — 1 and to the next to leading order in the derivative expansion. 
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1/m 

Figure 3: Critical exponents at various values of 1/m, including 1/m = 0. Empty 
and full black circles are respectively the 0(d°) and 0(d 2 ) values of v for N = 1. 
Empty and full black triangles are respectively the 0(d°) and 0(d 2 ) values of uj 
for N = 1. Finally the + and the x correspond respectively to the 0(d°) value 
of v and to the 0(d°) value of u for = 2. 
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1 2 



X 

Figure 4: The Wilson-Fisher fixed point derivatives of the potential (upper 
frame) and the eigenvectors corresponding to the exponent v (central frame) and 
those corresponding to the exponent oo (lower frame), for N — 1, 2, 3, 4, all to the 
0(<9°) order. 
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